REGULARITIES OF THE DISTRIBUTION OF ABSTRACT VAN DER 

CORPUT SEQUENCES 



WOLFGANG STEINER 

Abstract. Similarly to /3-adic van der Corput sequences, abstract van der Corput se- 
quences can be defined by abstract numeration systems. Under some assumptions, these 
sequences are low discrepancy sequences. The discrepancy function is computed explicitly, 
and the bounded remainder sets of the form [0, y) are characterized. 



1. Introduction 

Let {xn)n>o be a sequence with x„ G [0, 1) for all n > 0, and 

D{N, I) = i^{0 <n < N : x„ G /} - A^A(/) 

its discrepancy function (or local discrepancy) on the interval /, where A(/) is the length 
of /. Then, (x„)„>o is said to be a low discrepancy sequence if supj D{N, I) = O(logA^), 
where the supremum is taken over all intervals / C [0, 1). If D{N, I) is bounded in N, then 
/ is called a bounded remainder set. For details on discrepancy, we refer to |KN74l [DT97] . 
References to results on bounded remainder sets can be found in the introduction of |Ste06] . 

In |BG96l INin98al INin98bj . /3-adic van der Corput sequences are defined, and it is 
shown that they are low discrepancy sequences if /3 is a Pisot number with irreducible 
/3-polynomial. Recall that a Pisot number is an algebraic integer greater than 1 with all its 
conjugates lying in the interior of the unit disk. We refer to Section [3] for the definition of 
the /9-polynomial. In |Mor98t [lM04j . these results were extended to piecewise linear maps 
which generalize the /3-transformation. The proof in [Nin98al INin98bj relies on the fact 
that cylinder sets of the /3-transformation are bounded remainder sets if /3 is a Pisot number 
with irreducible /3-polynomial. Under the same conditions on /3, bounded remainder sets 
of the form [0,y), < y < 1, were completely characterized in [SteOGj : the /3-expansion 
of y is finite or its tail is the same as that of the expansion of 1. 

If /3 is a Pisot number, then the language of the /3-expansions is regular, which means 
that it is recognized by a finite automaton. Therefore, these /3-expansions are special cases 
of abstract numeration systems as defined in [LROll ILR02] , see Section [31 In Section [2|, 
we define van der Corput sequences related to more general abstract numeration systems. 
Theorem [T] in Section H] provides a new class of low discrepancy sequences. Finally, Theo- 
rem [2] in Section O characterizes the bounded remainder sets of the form [0,y) with respect 
to these abstract van der Corput sequences, generalizing the results in jSteOGj . 
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2. Definitions and first results 

Let {A, <) be a finite and totally ordered alphabet. Denote by A* the free monoid 
generated by A for the concatenation product, i.e., the set of finite words with letters in A. 
The length of a word w & A* is denoted by \w\. Extend the order on A to A* by the 
shortlex (or genealogical) order, that is to say v<wifv = woTv<w, where v < w 
means that either |f| < |w| or |f| = \w\ and there exist p,v',w' G A*, a,b & A such that 
V = pav', w = pbw' and a < b. 

According to |LR01] . the triple S = {L,A, <) is an abstract numeration system if L is 
an infinite regular language over A, and the numerical value of a word w G L is defined by 

ya\s{w) = E L : V < w}. 

If vals(w) = n, then we say that w is the representation of n and write rep_5(n) = w. 

Denote by A'^ the set of (right) infinite words with letters in A. It is ordered by the 
lexicographical order, that is to say t < u if t = u or t < u, where t < u means that there 
exist p G A* , a,b E A, t', u' G A"^ , such that t = pat', u = pbu' and a < b. 

Assume that the language L grows exponentially, with 

log#{t; G L : \v\<k} 
lim ; = logp > 0. 

Suppose that u G A^ is the limit of words w^''^ G L, i.e., every finite prefix of m is a prefix 
of w^''^ for all but a finite number of fc's. Then, the value of u is the real number 

/T\ ^ r val5(w('=)) 
1 valq(u) = lim -— — ■ — TTTTT, 

if this limit exists and does not depend on the choice of w^'^K Conditions assuring the 
existence of this value are given in |LR02j . see also Lemma [3] and its proof. Let 

= {u e A'^ : u= lim w^^'' for some w^*"'^ G L\. 
Since val^(u) G 1], we define the normalized value 

We extend this definition to finite words w E L which are prefixes of words in L^^ by setting 
(w) = (u), where u is the smallest word in L^^ with prefix w. Since we want to define a 
sequence without multiple occurrences of the same value, we set 

L' = {w E L : (w) 7^ (v) for every v E L with v < w}. 

Recall that the mirror image of a word w = W1W2 ■ ■ ■ Wk, Wj E A, is w = Wk ■ ■ ■ W2W1 
and that the mirror image of a language L is L = {w : w E L}. Now, we are ready to 
define the main object of this paper, abstract van der Corput sequences. 

Definition 1 (Abstract van der Corput sequence). Let S = {L,A,<) be an abstract 
numeration system, where L is a regular language of exponential growth, every word w E L 
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is the prefix of some infinite word m G L^, and the hmit in ([T]) exists for every u E L^. 
Then, the abstract van der Corput sequence corresponding to S is given by 

Xn = (w) with w = Tepg,{n), 
where S' is the abstract numeration system (L', A, < ). 

Thus, the set of values of an abstract van der Corput sequence is {x„ : n > 0} = 
{{w) : w & L} = {{w) : w G L'}, and the position of {w), w G L', in the sequence is 
determined by the shortlex order on the mirror image of L'. We need a number of further 
assumptions on the language L in order to get precise formulae for the discrepancy. All 
these assumptions are satisfied by the /3-adic van der Corput sequence when the language 
of the /3-expansions is regular, cf. Section [31 and by Example [T] at the end of this section. 

Let Al = {Q, A, T, go, F) be a (complete) deterministic finite automaton recognizing L, 
with set of states Q, transition function r : Q x A ^ Q, initial state go and set of final 
states F. The transition function is extended to words, r : Q x A* ^ Q, hy setting 
r(g, e) = q for the empty word e and r(g, wa) = T{T{q, w),a). A word w E A* is accepted 
by Al, and thus in L, if and only if r(go, w) G F. 

Definition 2 (Totally ordered automaton). A deterministic automaton {Q, A,T,qQ, F) is 
said to be a totally ordered automaton if there exists a total order on the set of states Q 
such that, for all q,r E Q, 

q < r implies r(g, a) < T{r, a) for every a E A. 

From now on, all automata will be totally ordered automata. Furthermore, the maximal 
state will be the initial state and every state except the minimal one will be final. (In case 
Q = F, where the automaton recognizes A*, we add a non-accessible state to Q.) W.l.o.g., 
the set of states will be Q = {0, 1, . . . , ci} for some positive integer d, and the order on Q 
will be the usual order on the integers, hence q^ = d and F = {1, . . . ,d}. Moreover, we 
will assume that r(0, a) = for every a E A, i.e., the state is a sink. 

Lemma 1. Let L <Z A* be recognized by a totally ordered automaton Al = [Q, A,T,d,Q\ 
{0}), with Q = {0, 1, . . . ,d} and r(0, a) = for every a E A. Then, L is recognized by the 
totally ordered automaton A^ = {Q, A, t,0,Q\ {0}), where 

(2) r(r, a) = G Q : r(g, a) + r > d} for every r E Q, a E A. 

In particular, we have r(0, a) = for every a E A. 

Proof. A deterministic automaton A' = {Q', A, r', qQ, F') recognizing L is obtained by 
determinizing the automaton which is given by inverting the transition function r, see e.g. 
|Sak09] . This means that 

• Q' is a subset of the power set V{Q), 

• T'{q', a) = {q E Q : T{q, a) G q'} for every q' E Q', a E A, 

• the set of final states in Al is the initial state of A', i.e., q'^ = {1, . . . ,d}, 

• the final states in A' are those elements of Q' which contain the initial state of Al, 
i.e., F' = {q' eQ' : dE q'}. 
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We show that Q' C {g^, g^, . . . , g^}, where q'^ = {r + 1, . . . ,d} (g^ being the empty set). 
Since Al is totally ordered and r(0, a) = 0, we obtain that 

t'{%, a) = {q E Q : T{q, a) > 0} = {r + 1, . . . , rf} = g^ for some r E Q. 

In the same way, we get, for every r E Q with g,'. E Q', that 

r'(g^, a) = {q E Q : T{q, a) > r} = {s + 1, . . . , d} = q'^ for some s E Q. 

This shows that Q' C {go,^!, . . . ,qd}- It is easy to see that A' is a totally ordered au- 
tomaton, with the order on Q' given by g.^ < g^ if g^ C g^, i.e., r > s. We clearly have 
F' = {gg, gi, . . . , n Q'. If we extend the set of states to {q^, q[, . . . , g^} (with possibly 
non-accessible states) and label the states hj d — r instead of g^, we obtain Ai- Therefore, 
Al is a totally ordered automaton, with Q ordered by the usual order on the integers. □ 

The next lemma provides a fundamental characterization of the words in a language L 
recognized by a totally ordered automaton Al = {Q, A, r, d, Q\{0}) with Q = {0,1, ■ ■ ■ ,d} 
and r(0, a) = for every a E A. 

Lemma 2. Let L, r, r be as in LemmaUl wi ■ ■ ■ Wk E A* , < j < k. We have Wi ■ ■ ■ E L 
if and only if T{d, Wi - ■ ■ Wj) + T{d, w^ - ■ ■ "U^j+i) > d. 

Proof. Let < j < k. With the notation of the proof of Lemma [H T{d,Wk- --Wj^i) = 
d — r can be written as r'^q'^, Wk - ■ ■ Wj+i) = q^ = {r + I, . . . ,d}. In Al, this means that 
Wj^i ■ ■ - Wk leads to a final state from the state g, i.e., r(g, Wj^i ■ ■ ■ Wk) > 0, if and only if 
q > r. Therefore, we have T{d, Wi - ■ ■ Wk) = r^r^d, Wi - ■ ■ Wj), Wj+i ■ ■ ■ Wk) > if and only if 
T{d, wi - ■ ■ Wj) > r, i.e., T{d, Wi - ■ ■ Wj) + T{d, w^ - ■ ■ Wj^i) > d. □ 

Remark 1. If we consider T{d, a) + ■ ■ ■ + r(l, a) as a partition of an integer, then r(d, a) + 
■■■ + r(l,a) is the conjugate partition, since T{d — r,a) = #{g E Q : T(q,a) > r}. 



E.g., if (r(d, a),..., r(l, a)) = (4, 2, 1, 0), then (r(ci, a),..., r(l, a)) = (3, 2, 1, 1), and the 



Next, we characterize the values of the abstract van der Corput sequence, under the 
assumption that the incidence matrix of the co-accessible part of Al is primitive. (A state 
g is co-accessible if T{q,w) E F for some w E A*.) 

Lemma 3. Let L be as in LemmaUl o,nd assume that Ml = E A : T{q, a) = ?"})i<5,r<d 
is a primitive matrix. Then, the normalized value exists for every u E and is given by 



corresponding Ferrers diagram is 



3 2 11 



4 
2 
1 




oo 



(3) 
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where /3 is the Perron- Frohenius eigenvalue of M^, [rji, . . . jTiaY is the corresponding right 
(column) eigenvector with rja = 1, tjq = 0, and u = U1U2 ■ ■ ■ with Uj G A for all j > 1. 

Proof. The definition and primitivity of tlie incidence matrix give 

(4) G A'' : r(g, t;) > 0} = (0, . . . , 1, . . . , 0)M^^(1, ...,!)* = cr^,/3^ + 0{p') 

witli constants c > and p < (3 sucli that every eigenvalue a 7^ /9 of Ml satisfies |a| < p. 
Due to the assumptions on Al, we have Ui---Uk G L for every k > 1. Similarly to 
[LR0n[LR02] . we split up 

{f G L : V <ui - ■■ Uk] = {f G L : |f I < A;} U [J [J {mi ■ ■ ■ Uj-iaw E L : w E A''^^. 

l<j<k a<Uj 

Since G L : \v\ < k} = Ej=o(c/3^ + ^(/^)) = + C(max(l, p)'=), we obtain 



valsiui ■■■Uk) 



1 /3 - 1 .^fma.x(l,p)''\ 



#{w G L : \v\ <k} 13 

Therefore, we have 

.u^. ^ valsiui ■ ■ ■ Uk) 1/3-1^ J 

valo u = hm —— — — = — H -; — > ej{u)p ■' , 

' j^{v e L : \v\ < k} (5 1^ jr[ 

and (m) = (/3vals(M) - - 1) yields (E]). □ 

As a last preparation for the study of the discrepancy of abstract van der Corput se- 
quences, we consider the language L' . 

Lemma 4. Let L he as in Lemma and assume that T{q, oq) > for every q > 0, where 
denotes the smallest letter of A. Then, L' consists exactly of those words in L which do 
not end with a^. 

Proof. We clearly have {t) < (u) if t < u, t,u E L^^. The primitivity of Ml implies that 
Tjq > for all g > 0. Therefore, (t) = (u) if and only if there exists no word u' G L^^ with 
t < u' < u. For V < w, v,w E L, this means that (v) = (w) if and only if f is a prefix of w 
and w is the smallest right extension of f in L of length \w\. Since r(g, ao) > for every 
g > 0, we have vQq E L for all k > 0. It follows that (v) = {w) with v < w, v,w E L, ii 
and only ii w = wa{j^'~''''. Thus, ti^ ^ L' if and only if w ends with oq. □ 

Example 1. Let Al = ({0, 1, 2, 3}, {qq, ai, 02}, r, 3, {1, 2, 3}) be the totally ordered automa- 
ton in Figure [Hon the left. The first words in L (in the shortlex order) are 



£, Qq, Cil, (3-2; (^QOiQi (^O^^li Oi0(^2i Oil^^Oi 0,lO'2i (^2O'0i 02C^2) 'j'-O'^O'^O) CLqCLqCLi., CLq(1qCL2, CLqCLiCLq, (lQ(lia,2, 
CiQ(l2CiQj CiQ(l2Ci2} CL^CiQCiQ^ (l\CiQCi\^ (l\CiQCi2} (^l(^2^0} (^1(^2^2} ^2^Q^0^ Ci2^Q^2^ Oj2^2^Q^ (^2^2^2' 
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ao,a2 




ao, ai,a2 



Figure 1. The totally ordered automata Al (left) and (right) of Example [T] 
The transition functions r, r and the incidence matrices Ml, Mj^ are given by 
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Recall that Ml = (#{a e A : r(g, a) = r})i<g,^<3, = (#{a G A : r(g,a) = r})i<,,^<3 

and that r can be calculated using Remark [H Thus, L is recognized by the totally ordered 
automaton A^ in Figured] on the right. The characteristic polynomial of Ml is — 2x^ — 
x + 1, the dominant eigenvalue is (3 ^ 2.247, and (?7i, ?72, ''73)* = (/3^ — 2/3, — /3^ + 3/3 — 1, 1)* ^ 
(0.555, 0.692, 1)* is a right eigenvector of Ml- Since — 2x^ — x + 1 is irreducible, it must 
be the characteristic polynomial of Mj^ as well. The conditions of Lemma H] are satisfied, 
hence the first elements of the abstract van der Corput sequence corresponding to L are 



Xo 

Xi = (0002) 



(e) = 0, xi = (ai) = — , X2 = (02) 



V3 + V2 



X3 = (aoOi) 



113 + V2 



X5 = (0102) = -5- + ^> ^6 = («2a2) 



?73 

+ ^?2 , m 



/3 



X7 = [aoaoai) = — , xs = (aiaoai) = — + — , xg = {aoaoa2) 



V3 + V2 
/33 ■ 



Xio = (010002) = ^ — ' ^11 = (02O0O2) 



^3 + ^2 ?73 

/3 



Xi4 



/ \ V3 V3 I \ V3 ^ V2 V2 

X12 = (000102) = ^ + 3:13 = (O0O2O2) = h 

/ .V3.V3.V2 I \ V3^V2 . ri2 ri2 

- = ^ + ^ + = (020202) = +p-2 + p-3- 



Xo 



Xg X3 



X12 X4, Xi3 Xi 



Xs Xio X5 Xi4 X2 



x\\ xe xi5 
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3. /3-ADIC VAN DER CORPUT SEQUENCES 

To obtain Ninomiya's /3-adic van der Corput sequences, consider a totally ordered au- 
tomaton Al = ({0, 1, . . . , d}, A, T, d,{l, . . . , d}) on the alphabet A = {0, 1, ... , B}, with 
integers bg E A, 1 < q < d, such that r(g, a) = d for all a < bg, and r(g, a) = if and only 
if a > 6g or g = 0. Assume that is primitive (which is the only interesting case since 
r(g,0) 7^ d implies bq = 0), and let /3 be its Perron- Frobenius eigenvalue. Then, we have 

= X] '7r(g,«i-«^_ia) = Vd = Uj for cvcry u = U1U2 ■ ■ ■ E L^, j > I. 

Let tit2 ■ • • be the maximal sequence in L^, i.e., tj = &r(d,ti - tj_i) for all j > 1. Since Al 
is a totally ordered automaton, T{d, ti - ■ ■ tj-i) < r((i, ti - ■ ■ implies tj < tk or tj = tk, 
T{d, ti ■ ■ ■ tj) < T{d, ti- ■ - tk), thus tjtj^i ■ ■ ■ < t^tk+i ■ ■ ■ , in particular tjtj+i ■ ■ ■ < tit2 ■ ■ ■ ■ 
By the maximality of tit2 ■ ■ ■ , we have Yl^i ^jP^'' ~ 1- The sequence tit2 ■ ■ ■ is called the 
infinite expansion of 1 in base /3 or quasi- greedy expansion of 1, cf. |Par60j . By the special 
structure of Al and the primitivity of Ml, we have {r((i, ti - ■ - tj) : < j < c?} = {1, . . . , d}, 
hence T{d,ti ■ ■ - td) = T{d,ti ■ ■ - tm) for some m < d. Therefore, ^1^2 • • • is eventually 
periodic with preperiod length m and period length d — m, which implies 



d y m \ 

7=1 ^ 7=1 ^ 



thus /3 is a root of the polynomial 

(5) {x^ - tix'^-^ tdX^) - (x™ - tix™-^ t^x°). 

If Al is the minimal deterministic automaton recognizing L or, equivalently, m and d — m 
are the minimal preperiod and period lengths of tit2 ■ ■ ■ , then (jS]) is called (3 -polynomial. 

If M = U1U2 ■ ■ ■ G L^^, then we have either u = tit2 ■ ■ ■ or mi ■ ■ ■ Uk-i = ti - ■ ■ tk-i, Uk < tk 
for some k > 1. Since T{d,ui ■ ■ - Uk) = d in the latter case, u E is equivalent with 
UjUj+i ■ ■ ■ < tit2 ■ ■ ■ for all j > 1, cf. [FarGOj . Therefore, u is either the greedy or the 
quasi-greedy /3-expansion of (u). Since wOO ■ ■ ■ is a greedy /3-expansion for every w E L, 
the abstract van der Corput sequence given by S* = (L, A, <) is exactly the /3-adic van der 
Corput sequence defined in [Nin98aj . Consequently, we call Al a (3 -automaton. 

Conversely, let ^1^2 ■ ■ ■ be the infinite expansion of 1 in base /3 > 1, i.e., the unique 
sequence of integers satisfying Yl^=i — ^ 00 ■ ■ ■ < tjtj+i ■ • ■ < ^1^2 ■ ■ ■ for all 

i > 1, cf. [ParGO] . Assume that tit2 ■ ■ ■ is eventually periodic. In particular, this holds 
when /3 is a Pisot number, see |Ber77l ISchSOj . Let d be the sum of the minimal preperiod 
and period lengths, and qj = #{1 < k < d : tktk+i ■ ■ < tj+itj+2 ■ ■ ■ } for j > 0. Then, the 
/3-automaton is given by bg^ = tj+i and T{qj,tj+i) = qj+i, see Example [2J 

Note that /3-adic van der Corput sequences were first considered in [BG96j for some 
cases in which L = L. In this case, the definition of the sequence is simpler beacause the 
/3-expansion of Xn is the mirror image of the expansion of n in a numeration system with 
respect to the linear recurrence corresponding to the /3-polynomial. In our notation, L = L 
is equivalent with t = t, see also |BY00l IKwo09] for a different characterization. 
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Example 2. Let /3 be the real root of — Ax^ — 2. Then, tit2 ■ ■ 
have d = 3, qo = 3, qi = 1, q2 = 2, q^ = 3, thus 63 = 4, 61 = 0, 62 
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= 401401 

1, and 



hence we 




This example is also given in Section 2.3 in |Ste06j . with different notation. The substitu- 
tion r in |Ste06j plays the role of the transition function r in this paper. 



4. Discrepancy function 

For a given abstract van der Corput sequence (x„)„>o and y G [0, 1], we study now the 
behavior of the function D{N, [0, y)) = #{0 < n < N : x„ < y} — Ny as — )■ 00. Since 
D{N, [y,z)) = D{N, [0,z)) — D{N, [0,y)), this determines the discrepancy of (x„)„>o. 

If L satisfies the conditions of Lemma [3l then every y G [0, 1] is the numerical value of 
some u & Li_j, see |LR02] . We have the following lemma, where 

Lemma 5. Let {xn)n>o be an abstract van der Corput sequence with L as in Lemma\^ 
For y = (u) with u = U1U2 ■ ■ ■ G L^, N > with rep^,(A^) = wi - ■ ■ Wi, we have 

e e 

j=l a<Uj k=j+l b<wi^ 

where L^j, = {f G : r(g, v) + r > d} = {v e : q + r(r, v) > d} and 

e 

C{N, n) = ^ #{6 < Wk : bwk+i ■ ■ ■ wia^ < UkUk+i ■■ ■ , ui- ■■ Uk-ihwk+i ■ --Wi E L}. 

k=l 

Since £ = 0{\ogN), we have C{N,u) = 0{hgN). 

Proof. We have to count the number of words v E L' with val^,('u) < A^, i.e., v < Wi - ■ - Wi, 
and {v) < y. As in the proof of Lemma HI we have wOq"'^' G L, thus {v) = (fOQ '^'). Since 
Wi > ao, V < Wi ■ ■ ■ Wi holds if and only if Oq"'^'?; < Wi ■ ■ ■ Wi. Therefore, we can count the 
number of words v E L (1 with v < Wi ■ ■ ■ Wi and (v) < y, instead of those in L'. 

The inequality (v) < y is equivalent with f Oq < u. Thus, we have to count the f in L 
of the form v = ui - ■ ■ Uj-iavj+i ■ ■ ■ Vk-ibwk+i ■ ■ -Wi with a < Uj, b < Wk, I < j < k < i, or 
V = U\ - ■ ■ Uk-ibwk+i ■ ■ - Wi with b < Wk, bwk+i ■ ■ ■ Wia^ < UkUk+i ■ ■ ■ , 1 < k < i. 

Lemma [2] yields that ui - ■ ■ Uj-iavj+i ■ ■ ■ Vk-ibwk+i ■ ■ ■ Wi E L ii and only if Vj+i ■ ■ ■ Vk-i G 
L7j~^ , h\, which provides the main part of the formula. The other words 

give C{N,u). We have i = 0{logN) since L grows exponentially and the same holds for 
L' by Lemma m Since A is finite, we obtain C{N,u) = 0{\ogN). □ 
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Similarly to |Nin98a[ INinQSbt ISteOGj , we assume now that the characteristic polynomial 
of Ml is irreducible. Let (32, ■ ■ ■ , (3d be the conjugates of the Perron- Frobenius eigenvalue 
(3i = (3. For any z e Q{(3), denote by ^ q(/5.) the image of z by the isomorphism 
mapping /3 to Similarly to (jlj), we have some constants 9i, . . . ,9d E Q{(3) such that 

d 

#Lj,, = (0,...,l,...,0)M|(0,...,0,l,...,l)* = 5^r/«^«/3f for 1 < g, r < 

i=l 

Note that (6'i, . . . , OdY is a right eigenvector of Mj^. Set Oq = and define 
(6) 1k{N) = E 07id,wr-wk+ib) for N = va\g,{w£ ■ ■ ■ wi), I < k < i, 

similarly to ej{u). Then, 



a<Uj b<Wf^ 1=1 



d e 



^ = E E = E E j:^f(^^Lr-.k,.^^t' = EE^i'(^)/5'"- 

fc=l fe<uife fc=l b<Wk *=1 *=1 fc=l 

For y = (m), we have thus 

e e oo 

D{N, [0, y)) = E E E E + ^(^' ^) - ^ E 

oo d / e £ \ 

= EE E f W7?(iv)/3r^'"' - E7i^^(^)/5'~'^^-(")/5"'' +^(^-«) 

j=l i=l \k=j+l k=l / 

oo / e d 

j=l V k=j+l i=2 

mm{j/) d \ 

- E Y.^l:\N)(3t'e,{u)(3-A+CiN,u). 

k=l i=l / 

Since the series converges absolutely, we can change the order of summation, and get 

e / k-l d 

D{N, [0,y)) = E EE7i^^(^)/5'''G?^(«)/5r^' - 

k=l \ j=l i=2 
oo d 



EE7?(^)/5'"'^^-w/5'M+^(^' 

j=k i=l 



U . 



The following theorem states that D{N, [0,y)) = 0{\ogN) if \(3i\ < 1 for 2 < z < d, i.e., 
if /3 is a Pisot number. The conditions on Ml are subsumed in the following definition. 
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Definition 3 (Pisot automaton). A deterministic automaton is said to be a Pisot automa- 
ton if the incidence matrix of its restriction to the states which are both accessible and 
co-accessible has one simple eigenvalue /3 > 1, and all other eigenvalues satisfy |a| < 1. 

Theorem 1. Let S = {L, A, <) be an abstract numeration system where L is recognized by 
a totally ordered Pisot automaton Al = ({0, 1, . . . , d}, A, r, d, {1, . . . , d}) with r(0, a) = 
for every a E A, r(g, oq) > for the minimal letter ao E A and every q > 0. Then, the 
corresponding abstract van der Corput sequence is a low discrepancy sequence. 

Proof. If Al is a Pisot automaton, then the characteristic polynomial of is irreducible, 
thus Ml is primitive. Therefore, the conditions of Lemma are satisfied, and D{N, [0,y)) 
is given by ([H]). Since 'jk{N) and ej{u) take only finitely many different values and i = 
0(\ogN), we obtain 

e 

D{N, [0,y)) = Y,0{1) + O(logiV) = O(logiV), 

k=l 

where the constants implied by the O-symbols do not depend on y. With D{N, [y,z)) = 
D{N, [0, z))-D{N, [0, y)), we get supj D{N, I) = 0(log N), and the theorem is proved. □ 

As Section |3]shows, the /3-adic van der Corput sequences defined by Pisot numbers /3 with 
irreducible /3-polynomial considered in |Nin98at INin98b] are special cases of the abstract 
van der Corput sequences in Theorem [H By Lemma [21 S = (L, A, <) is another, usually 
different, abstract numeration system satisfying the conditions of Theorem [T] whenever 
r(g, Oq) > for every g > 0, in particular when Al is a /3-automaton. The abstract van 
der Corput sequence considered in Example [T] is also a new low discrepancy sequence. 

5. Bounded remainder sets 
Under the conditions of Theorem [H ([7]) gives 

(9) DiN,[0,y)) 

oo / e d min(j,^) \ 

= E E E lk{N)e,{u)P''-^-']+C{N,u) + 0{l) 

j=l \ k=j+l i=2 k=l / 

for y = (u). If there exists some m > such that Uj = oq for all j > m, then ej{u) = 
for all i > m and C{N,u) is bounded. It follows that the interval [0, {v)) = [0, (fag)) is a 
bounded remainder set for every finite word v E L. 

For any /3-adic van der Corput sequence where /3 is a Pisot number with irreducible 
/3-polynomial, the bounded remainder sets [0,y) have been characterized in |Ste06j by the 
fact that the tail of the /3-expansion of y is either 0"^ or a suffix of tit2 ■ ■ ■ . With our notation 
and the /3-automaton defined in Section [3l this means that Yl'jLm+i = Vg some 

m > 0, g G In the more general case, we have the following partial characterization. 
Note that the proof is simpler than the one in [SteOGj . 
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Proposition 1. Let (x„)„>o be an abstract van der Corput sequence satisfying the condi- 
tions of TheoremUl u E L^. If there exists some m > 0, q E Q , such that 

where tq = tg,itg,2 ■ ■ ■ is the maximal sequence in with T{q, tq^i ■ ■ ■ tqj) > for all j > 1 
ifq>0,to = a^, eqj{tq) = T.a<t,,, VT{q,t,.i-tg,,_^a), thcn D{N, [0, {u))) is bounded. 

Proof. If g = 0, then ej{u) = for every j > m, and we have aheady seen that the 
boundedness of D{N, [0, (m))) follows from 

For g > 0, note that eq^2{tq)€-q,z{tq) ■■■ = eg',i(V)^9',2(V) ' ' ' ' where q' = r{q,tq^i), and 
that Ml is primitive. Therefore, we can assume that m is large enough such that has 
only positive entries. Hence, there exists some v E L (1 A"^ such that T{d, v) = q. Then, 
we have vtq G and 

em+l{vtq)em+2{vtq) ■■■ = €9,1 (tq)eg,2 (tg) " " " , 

thus D{N, [0, (u))) = D{N, [0, {vtq))) + 0(1) by ©. 

It only remains to show that D(^N, [0, {vtq))) = 0{1). If vtq = td, then {vtq) = 1 and 
D{N, [0, {vtq))) = 0. Otherwise, the successor v' of f in L has length m, and {vtq) = {v') 
since no sequence u' G L^^ satisfies vtq < u' < v'a^, hence 

D{N, [0, {u))) = D{N, [0, {vtq))) + 0(1) = D{N, [0, {v')) + 0(1) = 0{1). □ 

Clearly [y, z) is a bounded remainder set if both [0, y) and [0, z) have this property. For 
the converse, the following holds. 

Proposition 2. Let {Xn)n>a be an abstract van der Corput sequence satisfying the condi- 
tions of TheoremUl If D{N,I) is bounded, then X{I) G Q(/3). 

Proof. The proof is very similar to that of Theorem 1 in |Ste06] . Therefore, we only give the 
main steps. Define a substitution q i-> r(g, a^^o) • • ■ T^il, iq,mq), with {og^o, «g,i, • • • , iq,mq} = 
{a E A : T{q, a) > 0} and Uq^ < Qg^i < ■ ■ ■ < ag^m,, 1 < ? < c^, which plays the role of the 
substitution r in |Ste06j . Since T{d,ao) = d, we have d i-)- dw for some w G A*. Then, a 
continuous successor function on satisfying rep^,(n)aQ 1— >■ rep_g,(n + l)aQ is topologically 
conjugate to the successor function SonV defined in |Ste06j . see also |BR07j . If D{N, I) 
is bounded, then A(/) is an eigenvalue of the dynamical system (V, S), see Theorem 5.1 in 
|Sha78] and |Ste06] . By Proposition 5 in |FMN96j . these eigenvalues are in Q(/3). □ 

If y E Q(/3), then y = {u) for some eventually periodic sequence u, see |RS05j . Let p be 
the period length of ei{u)e2{u) ■ ■ ■ and m the preperiod length. From we get that 

e / d k-l 00 \ 

D{N, [0,y)) = Y,lY,^'^\N)Y,ef{u)f3t'~'-^,{N)Y,eAu)f3'-^^^^ 

k=l \ i=2 j=l j=k J 
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Set Uk = Yl^=k^A'^)P'^~^~^ ■ -^'-'^ k > m, we have 
ek{u)/3P'^ + --- + ek+p-i{u) 



Vk 



Y.^iy)Pt-''' = " ' , _ ,r + = + OW^) for 2 < z < rf, 

which gives 



D(iV, [0,y)) = CiN,u) - $^5^7?(iV)y? + ^(l). 

fc=i 't=i 

Recall the definition of 7fc(A^) in ([6]), and set 

d 

(10) Cr-(^) = ^r^^^^ ioTreQ,ze Q(/3). 

i=l 

Then, we have 

e 

(11) D{N, [0,y)) = C{N,u) - Cnd^^r-^.^AVk) + 0{l). 

k=l b<Wk 

Note that 







1 if g + r > (i, 



(12) = - ^ d,,. 

The main result of this section is the following theorem. 

Theorem 2. Let (x„)„>o he an abstract van der Corput sequence defined by an abstract 
numeration system S = {L, A, <), where L is recognized by a totally ordered Pisot automa- 
ton Al = ({0, 1, . . . , d}, A, T, d,{l,..., d}) with r(0, a) = for every a & A, T{q, Oq) > q 
for the minimal letter G A and every q E {1, . . . , d — 1}. Then, D{N, [0, (u))), u G L^, 
is bounded in N if and only if {u) G Q(/3) and there exists some m > such that 

1 if vQq < UkUk+i ■ ■ ■ and Ui - ■ ■ Uk-iV G L, 
else, 



(13) CT{d,v){yk) 



for every v G A*, k > m, where yk = '^JLk^ji'^)^ ^ ^, ^ "is as in ( as in ^^). 

Proof. We show first that the conditions on -u G are sufficient for the boundedness of 
D{N, [0, (u))). Note that the assumption r(g, ao) > g for < g < d is not used here, but 
only the weaker assumption r(g, ao) > of Lemma HI For = val^, (w£ ■ ■ ■ Wi), we have 

k=l b<Wk k=m+l ^ ■' 

= C(iV,w) + 0(l), 
thus ^(A^, [0, (m))) = C(l) by ([11]). 
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Now, suppose that D{N, [0, (u))) = 0{1) for u G L^. By Proposition |2l we have 
{u) G Q(/3), thus u is eventually periodic by |RS05j . Let m > 0, p > 1 be such that 

U = Ui - ■■ UmiUm+l ' ' ' Um+pT and r(d. Ml ■ ■ ■ Um) = T{d, Ml ■ ■ ■ Um+p) ■ 

Consider v E A* and k > m. li v ^ L, then ( fT3l) holds since Cr(d,v){yk) = CoiVk) = 0. 
If f = v'ao, then ( !T3|) holds for v if and only if it holds for v'. Therefore, we can assume 
V G L'. If V is not the empty word, then we define integers 

Nh = val^, {{val)VQ-^) for h>0, 

with g > max{d — l,p) such that g + |f | is a multiple of p. Set v = Vk ■ ■ ■ Vi, with k < i. 
We show that 

(14) C(A^;„ n) = (/i + l)C(iVo, ^^) for all h>0. 

If Mm+i ■ ■ ■ Um+p = flo, then C{Nh, u) = and f lT^ holds. If Um+i ■ ■ ■ Um+p > cto' then g > p 
implies that bvj+i ■ ■ ■ f^ag < UjUj+i ■ ■ ■ , k < j < i, if and only if bvj^i ■ ■ ■ ve{aQv)^aQ < 
UjUj+i ■ ■ ■ . Since Al is totally ordered, the assumption r(g, ao) > q for < q < d implies 
T{d, Oq) = d and r(g, Qq) = c? for all g > 0, j > c? — 1, in particular r(g, Oq) = (i. Therefore, 
we also have Mi ■ ■ ■ Uj^ibvj+i ■ ■ ■ V£ E L ii and only if mi ■ ■ • Uj^ibvj+i ■ ■ ■ vi{aQv)^ G L. Since 
(7 + |f I is a multiple of the period length of u, f|T^ follows. 

Since r(g, Oq) = d for all g > 0, we also have r(g, Oq) = for all g > 0, and thus 
C7{d,{v4)''^e-vj+ib){yj) = C,T{d,vr-v,+^b){yj)- By ([II]), we get 

D(iV,, [0, {u))) = {h + l) {c{N,, u)-Y,Y. Cr(d,.,....,+,6)(y,)) + 0{l). 

^ j=k b<Vj ' 

Therefore, D[N , [0, (m))) = C(l) implies that 

•^■^^^ 5Z 5Z Cf(d,t;r-»^j+i6)(2/i) = C'(^o, m) for A^o = val_5,(wf ■ ■ ■ Vka\~^), k > m. 

j=k b<Vj 

We will show that f|T3|) holds, by considering f|l5l) for the integer A^q defined by the 
successor of v in L'. Let first v be the empty word, and A^^q = val^,(t>^aQ~^) with such 
that T{d, b) = for all 6 G A with < b < v'^. Then, (ITB]) gives 



CdiVk) = Cr(d,ao)(2/fc) = <^(^o,m) 



1 if < UkUk+i ■ ■ ■ , 
else, 



thus (fl3|l holds if f is the empty word. Assume next that there exists some v'^ > Vk 
with T{d,V£ ■ ■ -ffc+if^) > 0. Choose v[ such that T{d,ve ■ ■ -Vk+ib) = for all 6 G A with 
Vk < b < v'f,, and set A'q = val^,(f£ ■ ■ ■ Vk+iv'f^aQ~^). Then, f|T3l) holds because of 

CT{d,VfVk){yk) = Y, (T{d,ve-Vk+ib){yk) - ^ C7{d,ve-Vk+ib){yk) 

= C{N'^,u)-C{N^,u 



b<v'^ b<Vk 



1 if Vk - ■ ■ VittQ < UkUk+1 ■ ■■ , ui - ■■ Uk-iVk ■■ - ve e L, 
else. 
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Now, we proceed by induction. We know that f|T3l) holds if v is the empty word. Assume 
that f[T5]) holds for all words v of length i — k, and consider 

d 



C7{d,ve-v^:+^)(yk+l) - ^^^{d,vi-vk+r){f^yk ~ ^kiu)) 



i=l 
d 



i=l i=l b<ui^ 

Let Vk be the maximal letter with r(c?, vg - ■ ■ Vk+iVk) > 0. Since (^i, . . . , OdY is an eigenvector 
of Mj^ and we already know that (T(d,vf-vk+ib){yk), b < v^, is given by f|T3|) . we have 

d d 

f^i^T{d,ve-Vk+i)yk = '^2'^2^T{d,ve-Vk+ib)yk = CT{d,VfVh+ib){yk) 
i=l i=l beA feeA 

= Cr(d,i)r-'^fc)(^fe) + < ^fc • ■ ■ ■ ^^'^0 < UkUk+1 ■■■ , ui - ■■ Uk-ihvk+1 ■■■v^e L}. 

Using this equation, 

d 

Ed}) \^ {i) = \^ 4ir^ 

^T(d,VfVk+i) / J iT{d,ux-Ui,_xb) / J Tr'^T{d,ui---Uk-ib),T{d,vi---Vk+i) 

1=1 fe<Mfc b<u^: 

= #{6 < Mfe : ui ■ ■ ■ Uk-ibvk+i ■■■Vie L} 

and the induction hypothesis, f|T6|) yields 

Cr(d,i;r-^'fe)(2/fc) = Cf -^-fc+i) (^/fe+l) + < ^fe : Ui • ■ ■ ■ • • t;^ G L} 

- #{6 < t;fc : bvk+i ■ ■ ■ VittQ < UkUk+i • ■ ■ , Mi ■ • • Uk-ihvk+i ■■■Vi>eL] 

1 if Vk+i ■ ■ ■ Vfa'^ < Uk+iUk+2 ■ ■ ■ and Mi ■ ■ ■ UkVk+i ■ ■■ve e L, 

else, 

{1 if < Uk and Mi ■ ■ ■ Mfc-it;^ ■ ■ ■ G L, 
-1 if > Uk, Vk+i ■ ■ ■ vettQ < Uk+iUk+2 ■■■ , ^i ■ ■ ■ UkVk+i ■■■v^e L, 
else, 

1 if f fc ■ ■ ■ f^flo < UkUk+i ■ ■ ■ and Mi ■ ■ ■ Uk-iVk ■ ■■viE L, 
else, 

where we have used that Mi • ■ ■ Uk-ibvk+i ■ ■ ■Vi can be in L only if 6 < f a;. Therefore, (ITS]) 
holds for all words f of length £ — /c + 1, and Theorem [2] is proved. □ 

In the case of /3-adic van der Corput sequences, all bounded remainder sets satisfy the 
conditions of Proposition [T], see [SteOGj . The following example shows that this is probably 
not true in the more general setting, i.e., that there might be sequences u satisfying the 
conditions of Theorem [H but not those of Proposition [1] However, Al is not a Pisot 
automaton in this example, and we have not found an example where Al satisfies the 
conditions of Theorem [H 
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Example 3. Let Al = ({0, 1, 2, 3, 4}, {a^, ai, 02, 03}, r, 4, {1, 2, 3, 4}) be the totally ordered 
automaton with r, r given by the transition tables 



T 


ao 


ai 


^2 


as 


T 


ao 


ai 


(32 


as 
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2 


1 
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1 
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2 
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2 
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4 
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3 


1 


4 


4 


4 


4 


2 


4 


4 


4 


4 


3 



If u = a^aQa"^, then r(4, Ui ■ ■ ■ Uk-i) = 3 and thus €k{u) = 27/4 for all A; > 3. We obtain 
yk = Vi-m + Vi, which implies UiVk) = C,2{yk) = 1, CsiVk) = CiiVk) = by ([12]). It can be 
easily verified that f lT3|) holds for all f G A*, /c > 3, but the conditions on u of Proposition [1] 
are not satisfied. 

We conclude by the remark that the boundedness of D{N,I) is not invariant under 
translation of the interval, i.e., D{N, [y,z)) can be unbounded if [0,z — ?/) is a bounded 
remainder set and vice versa, see jSteOGj . 
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